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Strain has been widely used to manipulate the properties of various kinds of materials, such as ferroelectrics, 
semiconductors, superconductors, magnetic materials, and "strain engineering" has become a very active 
field. For strain-based information storage, the non-volatile strain is very useful and highly desired. 
However, in most cases, the strain induced by converse piezoelectric effect is volatile. In this work, we report 
a non- volatile strain in the (001) -oriented Pb(Mgi/3Nb2/3)03-PbTi03 single crystals and demonstrate an 
approach to measure the non-volatile strain. A bipolar loop-like S-E curve is revealed and a mechanism 
involving 109° ferroelastic domain switching is proposed. The non-volatile high and low strain states should 
be significant for applications in information storage. 



Strain has been widely used to manipulate the properties of various kinds of materials^"^, such as ferro- 
electrics\ semiconductors^, superconductors^, magnetic materials^, etc., and "strain engineering" has 
become a very active field\ One popular way to introduce strain is via lattice misfit between the epitaxial 
film and the substrate, and dramatic variation of strain can be realized by selecting suitable film and substrate\ 
However, in this case, the strain can not be tuned continuously and reversibly, and different substrates may induce 
other effects besides the strain effect when comparing the properties of films grown on different substrates. In this 
regard, converse piezoelectric effect is very useful for introducing continuous and reversible strain by electric 
fields. Recently the strain -mediated multiferroic heterostructures, composed of magnetic thin films and ferro- 
electrics (FE), have attracted much attention due to the magnetoelectric coupling between magnetic thin films 
and FE^"^^. In these multiferroic heterostructures, strain is produced via the converse piezoelectric effect of FE and 
can be modulated continuously and reversibly. This kind of strain becomes a primary pathway to control the 
electronic transport and magnetic properties of the magnetic films in the multiferroic heterostructures, and 
"straintronics" has been proposed^^ Owing to the superior piezoelectric properties'^, relaxor ferroelectrics such 
as (l-x)Pb(Mgi/3Nb2/3)03-xPbTi03 (PMN-PT) single crystals, which are binary solid solution system (x repre- 
sents the nominal composition of PbTi03), have been widely used in the strain-mediated multiferroic hetero- 
structures, besides their applications in ultrasonic devices, telecommunication, medical imaging'^. So far, most 
reported in-plane strain- electric field (S-E) behavior of the (001) -oriented PMN-PT is a symmetric butterfly-like 
curve without remanent strain'^, i.e. the electric-field induced strain is volatile. However, non-volatile strain is 
more important and highly desired for the strain-based information storage in the multiferroic heterostructures'^. 
For electric-field induced reversible non-volatile strain, there are only a few reports on FE materials with remnant 
strain states in the S-E curves, which can be classified into two types. For type I, the polarization -electric field (P- 
E) hysteresis loop is offset about the origin as shown in polycrystalline PbZri.xTix03 due to an internal field which 
is generated by different kinds of defects^^. Because the correlation between strain and polarization, the S-E curves 
also deviate from the symmetric butterfly-like curves and show an asymmetric behavior. In this case, there is 
indeed remanent strain at zero electric field, but the asymmetric behavior of the butterfly- like curves are deter- 
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Figure 1 | Experimental configuration, ferroelectric property and two different methods for measuring strain, (a) Schematic of the sample and the 
experimental configuration, (b) P-E loop of PMN-PT. (c) and (d) Schematic of the continuous and pulsed measurement methods for S-E curves, 
respectively. The strain measurement starts at the end of interval as pointed by the solid blue triangle symbols. 



mined by the polarity of the space charges^ \ accompanied with 
fatigue^^'^^. For type II, Wu et aU"^'^^ recently obtained the remanent 
strain in the (Oil) -oriented PMN-PT single crystals. When a critical 
reversed electric field was applied, the polarizations reoriented to the 
four possible in-plane <111> directions, which acted as a meta- 
stable state during the polarization reversal process. After decreasing 
the electric field to zero, hysteresis occurred for both P-E and S-E 
curves and two different strain states arise. However, a critical 
reversed electric field is required to induce the effect in these systems. 
If the electric fields applied to PMN-PT are larger or smaller than the 
critical field, the remanent strain will decrease dramatically^^. Since 
the critical reversed electrical field is a bit smaller than the coercive 
field, the cycle of the electric field applied to the (01 l)-oriented PMN- 
PT is actually unipolar as mentioned in Wu's paper^^, while for the 
bipolar case, there is no remanent strain^'*. In this paper, we report a 
non-volatile strain in the (001) -oriented PMN-PT single crystals. A 
bipolar loop-like S-E curve is revealed and a mechanism involving 
109° ferroelastic domain switching is proposed. We also put forward 
and demonstrate an approach for measuring the non-volatile strain. 

Results 

Methods for strain measurement and loop-like strain behaviors in 
PMN-PT single crystals. Figure la is the schematic of the sample 
and the experimental configuration. (001) -oriented (l-x)PMN-xPT 
(x = 0.28, PMN-PT with 28% of PT) single crystals were used in our 
work. The single crystal orientation and microstructure of the 
polished surface were characterized as shown in Fig. SI in 
Supplementary Information. Figure lb shows the P-E hysteresis 
loop for PMN-PT measured by using a Radiant Technologies 
Precision Premier II system. The in-plane S-E curves were 
measured by using the continuous method as shown in Figure Ic, 
which has been commonly used in the literature^^'^^. We measured 
the in-plane strains along the [110] and [1-10] directions of the (001) 
plane, whose converse piezoelectric responses are more sensitive^^. 
The S-E curve along the [110] direction is shown in Figure 2a. In 
contrast to the normal symmetric butterfly-like curve^^, it shows an 
asymmetric behavior with two remanent strain states at zero electric 
field. The in-plane S-E curve along the perpendicular direction 
[1-10] also shows an asymmetric butterfly-like behavior (Figure S2 



in Supplementary Information) and is complementary with that of 
the [110] direction, providing a clue to the origin of the phenomenon 
as discussed later. 

As shown in Fig. lb, the hysteresis loop of our sample is symmetric 
about the origin, which means that the internal field induced by 
defects is negligible. This is different from the aforementioned type 
I asymmetric S-E curve, which is related to the asymmetry of the P-E 
hysteresis loop. Moreover, unlike the unipolar S-E curves for the 
aforementioned type IP^, the S-E curve of our sample is bipolar. In 
addition, the electric fields apphed in our work (±8 kV/cm) are 
much higher than the coercive field (1.2 kV/cm) characterized in 
the P-E loop (Fig. lb). So the sample is in the saturated polarization 
states, which is different from the aforementioned type II case. Thus, 
a different mechanism for the asymmetry of S-E curves should be 
involved in our samples. 

In order to get the non -volatile part in the asymmetric butterfly- 
like curves, we proposed an approach as shown in Fig. Id, which is 
denoted as the pulsed measurements method. The detailed descrip- 
tions for the continuous and pulsed measurements are in the experi- 
mental section. In both measurements, we begin to measure the 
strain at the end of the interval. The difference is that electric voltage 
is always applied to the sample in the continuous measurement, 
while no electric voltage during the interval for the pulsed measure- 
ment. Thereby, in the pulsed measurement, the strain is measured at 
zero electric filed and reflects the remanent strain. Figure 2b shows 
the S-E curve obtained by the pulsed measurement. Interestingly, it 
exhibits a bipolar loop -like behavior, similar to the P-E hysteresis 
loop (Fig. lb). The loop-like S-E curve suggests that one can easily 
control the strain states by electric fields. For the perpendicular 
[1-10] direction, we also got the asymmetric butterfly- like and 
loop-like S-E curves (Fig. S2), which exhibit a reverse behavior com- 
pared to those of the [110] direction. Figure 2c shows the S-E curve 
obtained by subtracting the strain values in the loop-like curve 
(Fig. 2b) from those in the asymmetric butterfly-like curve 
(Fig. 2a) at the same voltage point by point. The S-E curve in 
Fig. 2c is a symmetric butterfly-like curve, which is common in the 
previous reports^^. Similar results were also obtained along the [1-10] 
direction (Fig. S2). Based on the above results and analyses, it can be 
deduced that the observed asymmetric butterfly-like S-E curves are 
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Figure 2 | Electric-field tuned strain behavior along [110] direction, (a) and (b) S-E curves along the [110] direction measured by continuous and 
pulsed measurements, respectively, (c) The curve obtained by subtracting the strain values in the loop-like S-E curves (Figure 2b) from those in 
asymmetric butterfly curves (Figure 2a). (d) The stable and remarkable high/low remanent strain states achieved by switching the polarity of the pulsed 
electric field. 



composed of the symmetric butterfly-like and the loop-like curves. 
Stable and remarkable high/low remanent strain states can be 
achieved by switching the polarities of the pulsed electric fields as 
shown in Fig. 2d for the [110] direction and Fig. S2d for the [1-10] 
direction. It can be seen that the pulsed electric field polarity depend- 
ence of the high/low strain states for these two directions is converse, 
which is consistent with the complementary behavior in their S-E 
curves (Fig. 2 and Fig. S2). Instead of a specific value or range of 
electric field required in the type II case, any fields larger than the 
coercive field of the sample can induce the remarkable high/low non- 
volatile strain states. Moreover, as the sample is in the saturated 
polarization states, the non-volatile strain states were stable and 
didn't show an obvious relaxation with time, which is in contrast 
to the behavior of the (Oll)-oriented PMN-PT^^ And the high/low 
strain states can be switched by pulsed electric fields for thousands of 
times without obvious fatigue effect (Fig. S3 in Supplementary 
Information). 

Proportions of ferroelastic switching in PMN-PT single crystals 
with different S-E behaviors. It should be mentioned that not all the 
samples exhibit the asymmetric butterfly-like S-E curves. We also 
obtained the symmetric butterfly-like curve as shown in Fig. S4a in 
Supplementary Information, similar to that reported by Thiele^^, in 
samples with the same nominal composition as those showing the 
asymmetric butterfly-like behavior. Figure S4b shows the S-E curves 
measured by the pulsed measurement. Obviously, there is no 
remanent strain at zero field in this kind of samples. Hereafter, we 
denote the samples with the asymmetric and symmetric S-E curves as 
sample A and sample S, respectively. To identify the differences 
between these two kinds of samples, we carried out XRD-RSM 
with in situ electric fields since it can provide straightforward 
information on the microstructure of the samples and the 
switching category of FE domains. PMN-PT samples are in the 
rhombohedral phase with eight spontaneous polarizations along 
the <111> directions, denoted as rl^/rl", r2^/r2", r3^/r3~, 
r4^/r4" as shown in Fig. 3a and 3b. When the polarity of electric 
field was changed, the spontaneous polarizations can be switched by 
71°, 109° and 180°, respectively. Polarization variants changing from 



rl^/r3^ to r2"/r4" or r2^/r4^ to rl"/r3" (vice versa) correspond to 
the 109° switching, while polarization variants changing from rl^ to 
r3" or r3^ to rl" (the same case for r2/r4) correspond to the 71° 
switching, and polarization variants without any category change 
(e.g. rl^ to rl" and so on) correspond to the 180° switching. 
Figure 3c and 3d show the reflections around the (113) peak for 
sample A under —8 kV/cm and +8 kV/cm, respectively. The 
distortions corresponding to polarization variants r2 and r4 are 
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Figure 3 | The schematic diagram of the polarization variants and XRD- 
RSMs around (113) reflections, (a) and (b) Configurations of the (113) 
crystal face (translucent light blue plane) and polarization variants of 
rhombohedral phase under negative and positive electric fields, 
respectively, (c) and (d) (113) reflection peak of XRD-RSMs for sample A 
with in situ electric fields of — 8 kV/cm and + 8 kV/cm, respectively. 
Ellipses with different colors in RSMs outUne the contours of the deduced 
spots. 
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Table 1 | Proportions of different polarization variants under pos- 
itive and negative electric fields, respectively 

E(kV/cm) rl r2/r4 r3 Correlation Coefficient 

-8 41.7% 55.5% 2.8% 0.96 

+ 8 17.1% 35.6% 47.3% 0.94 



degenerate in the (113) reflection peak, so there are three spots 
corresponding to rl, r3, and r2/r4, respectively. According to the 
lattice parameters in the rhombohedral structure^^ and the 
analytical method^°, we analyzed the positions and intensities of 
the spots (see SI for details), which can give the percentages of 
different polarization variants under positive and negative electric 
fields, respectively. Table 1 shows the results for sample A, suggesting 
that the proportion of r2/r4 under +8 kV/cm is about 20% less than 
that under —8 kV/cm. Since the change between r2/r4 and rl/r3 
from positively poled cases to negatively poled cases (vice versa) 
corresponds to the 109° switching, the 20% reduction of r2/r4 
proportion represents the net proportion of 109° switching. For 
comparison, we also measured the reflections around the (113) 
peak with in situ electric fields for sample S. The results are shown 
in Figure S4c and S4d. Similar to sample A, the (113) reflection peak 
also shows three spots. We analyzed the data using the same method 
as for sample A and the results are shown in Table SI. The 
percentages of different polarization variants are almost 
unchanged with the polarities of electric field. This means that 
there are 71/180° switching, but almost no net 109° switching 
since the percentage of r2/r4 is nearly unchanged under positive 
and negative electric fields. 

The crucial role of the net 109° switching for the non-volatile 
strain. In sample A, there are not only 71/180° switching but also 
about 20% of net 109° switching. While in sample S, there are only 
71/180° switching and almost no net 109° switching. Combined with 
the complementary property in the loop -like curves in sample A 
(Fig. 2b and S2b) and the symmetric butterfly-like curves in 
sample S (Fig. S4a), it can be concluded that the loop-like S-E 
curve is related to the net 109° switching. In contrast, without net 



109° switching, the S-E curve exhibits the butterfly-like behavior. 
Figure 4 is the schematic diagram of the 109° switching, 71/180° 
switching and corresponding changes of distortion/strain. For 
polarization switching from to P~, we take rl^ and r3^ as an 
example to show the 109° switching (right) and 71/180° switching 
(left) as shown in the upper panel of Fig. 4. It can be seen that the 
projections of polarizations to the (001) plane before and after the 
109° switching are perpendicular, while they are parallel before and 
after the 71/180° switching. As a result, the in-plane tensile 
distortions related to the polarizations are perpendicular before 
and after the 109° switching and parallel before and after the 
71/180° switching (lower panel of Fig. 4). It is also noted that 
change of the tensile distortions along the [110] and [1-10] 
directions are complementary for the 109° switching. This is the 
origin of the complementary loop -like curves of sample A along 
the two perpendicular directions. In sample S, there are only 71/ 
180° switching and the S-E curves show a symmetric butterfly- like 
behavior without remanent strain. 

Discussion 

An interesting question arises as why sample A and sample S, with 
the same nominal composition, show such a big difference in the S-E 
curves, related to the presence and absence of net 109° switching in 
sample A and sample S, respectively. It is likely that defects/defect 
configurations in the crystals may play an important role. As 
reported by Fu et al.^\ polarization rotation along the different paths 
needs to overcome different energy barriers. And it has been shown 
that defects can strongly affect the barrier for polarization rotation, e. 
g. defects make the 180° polarization flip become a feasible process^^. 
It was also shown that defect configurations play an important role in 
the polarization switching^l These may shed light on our case as 
follows. Some defects/defect configurations in sample A lowers the 
energy barrier between the two polarization configurations with the 
polar axis rotating by 109°, which is crucial for the loop-like S-E 
curves. In contrast, it is not the case for sample S and 109° domain 
switching is not favored due to the required high energy. The pres- 
ence of key defects/defect configurations maybe due to details during 
the crystal growth process. Further study on this issue is deserved. 




SCIENTIFIC REPORTS | 4:4591 | DOI: 1 0.1 038/srep04591 



4 



In conclusion, we report a non-volatile strain in the (001) -oriented 
PMN-PT single crystals. A bipolar loop-like S-E curve is revealed and 
a mechanism involving 109° ferroelastic domain switching is 
proposed. We also put forward and demonstrate an approach for 
measuring the non -volatile strain. As a comparison, the non- volatile 
strain in the (01 l)-oriented PMN-PT reported previously is depend- 
ent on a critical reversed electric field and decreases dramatically if 
the electric field deviates from the critical value. While the non- 
volatile strain in the (001) -oriented PMN-PT in our work can be 
obtained by any fields larger than the coercive field and maintains 
a stable state. The non-volatile high and low strain states should be 
significant for applications in information storage. 

Methods 

Our samples are one- side-polished (001) -oriented PMN-PT with the edges cut along 
the pseudocubic [110] and [1-10] directions, and dimensions of 5 mm X 5 mm X 

0. 5 mm. Au layers of 100 nm thick were deposited on both sides of the samples as 
electrodes. The strain gauge was pasted on the sample surface with glue M-Bond 610, 
and kept at 120°C for 2 h to strengthen the paste effect. 

Two kinds of measurements for S-E curves are shown in Fig. Ic and Id, which are 
named as continuous and pulsed measurements, respectively. In the continuous 
measurement, the voltage applied to the sample is increased with a step of 20 V. After 
a period of delay, an interval is followed. At the end of the interval, the measurement 
begins, and then the voltage is changed for the next cycle. So the voltage is always 
applied. In the pulsed measurement, the voltage is also increased with a step of 20 V. 
After a period of delay, the voltage is removed for an interval. At the end of the 
interval, the strain measurement starts. After that the voltage is changed and the 
measurement goes into the next cycle. We carried out the XRD-RSM experiments 
with in situ electric fields of ± 8 kV/cm using the synchrotron radiation x-ray source. 
All the measurements were carried out at room temperature. 
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